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Abstract: The lowest-lying glueball masses are computed in SU(JV) gauge theory on a spacetime 
lattice for constant value of the lattice spacing a and for N ranging from 3 to 8. The lattice spacing 
is fixed using the deconfinement temperature at temporal extension of the lattice Nt = 6. The 
calculation is conducted employing in each channel a variational ansatz performed on a large basis 
of operators that includes also torelon and (for the lightest states) scattering trial functions. This 
basis is constructed using an automatic algorithm that allows us to build operators of any size 
and shape in any irreducible representation of the cubic group. A good signal is extracted for the 
ground state and the first excitation in several symmetry channels. It is shown that all the observed 
states are well described by their large N values, with modest 0(1/N 2 ) corrections. In addition 
spurious states are identified that couple to torelon and scattering operators. As a byproduct of 
our calculation, the critical couplings for the deconfinement phase transition for N = 5 and N = 7 
and temporal extension of the lattice Nt = 6 are determined. 
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1. Introduction 

Non-Abelian SU(iV) gauge theories with fcrmions in the fundamental representation have a simpler 
diagrammatic expansion in the limit in which the number of colours N goes to infinity with the 
product A = g 2 N held fixed pL 0, E| (see e.g. 0, p| for reviews of the original ideas and some more 
recent developments). Due to their simpler structure that can nevertheless shed invaluable insights 
on non-perturbative phenomena in non-Abelian gauge theories, SU(iV) gauge theories in the large 
N limit play a central role in the gauge-gravity correspondence and have become the subject of a 
line of numerical investigations on the lattice. In addition to determining values for observables in 
the large N limit, lattice calculations provide their corrections at finite N, which can be expressed 
as a power series in 1/N for the quenched theory and in 1/N in the dynamical case; the emerging 
picture is that, at least for the quenched theory, only the leading correction of 0(1/N 2 ) is sufficient 
to describe the system at any finite value of N bigger than two at a level of accuracy of the order 
of a few percents. These calculations clarify quantitatively the meaning of the statement that the 
physical case N = 3 is close to N — oo. Numerical calculations at large N have been recently 
reviewed in p, R, pi , with a review of earlier numerical works provided in H . 

If the lattice legitimates the use of the large N limit to describe the physics for N = 3, 
analytical approaches based on large N ideas are of great help to model the lattice data. In order 
to assess the reliability of various analytical methods based on the large N framework (which often 
have to resort to other approximations in addition to taking the large N limit), it is important to 
compare their predictions to the lattice data for observables that are well under control in both 
approaches. The glueball spectrum in the pure Yang-Mills theory is one of the easiest observables to 
compare. Analytical results for glueball spectra in the large N limit are available from calculations 
based on different approaches, in particular in various backgrounds in the gauge/string duality 
framework (l0| [yj [l2|, [13| [lj, [jjj [u| [17|, Q , with a variational ansatz for the groundstate of the 



Yang-Mills theory |T3, EG, Ell and in the light-cone framework [E2l E3l 24 . On the lattice side 



glueballs in the SU(2) and mostly SU(3) Yang-Mills theory have been the subject of extended 



studies for quite a long time |25j, Efl, E7K LM [29, 30 , 31j , and recent investigations have provided us 



with reliable determinations of the spectrum in all the J PC channel s |32|, p3|. Reliable numerical 
results at larger N have become available only more recently |34|, ^5|, |(| f37| |8| |3^, 0, 0, ff2| an( i 
mostly concern the groundstate and the first excitation in the ++ and 2 ++ channels. 



In general, hadronic states on the lattice are computed with a variational ansatz, and in order 
to reliably extract the energy of excitations the variational basis must be large. Building a large 
variational basis in all the symmetry channels, corresponding to the quantum numbers of the 
hadrons for which we want to extract the mass spectrum, proves to be quite complicated from 
a technical point of view (see e.g. Refs. ^, [J4J for a discussion in the context of the meson 
spectrum). For this reason, so far the bulk of the effort for extracting reliable single-particle 
spectra in the glueball sector has been put on the physical case N = 3, with the calculations at 
larger N focusing mostly on much fewer states. Given the theoretical importance of determining 
accurately the behaviour of the theory at large N, it would be useful to have a more comprehensive 
calculation also in that case. 

Our work aims to fill this gap. In this paper we will provide the first determination of the 
large N glueball spectrum (obtained with an extrapolation including values of N up to eight) in 
several irreducible representations of the lattice rotational group and for both values of parity and 
charge conjugation. We find that, also for the states studied for the first time in this work, a 
small 1/N 2 correction to the N = oo limit describes the data with an accuracy of a few percents 
all the way down to N = 3. Progress over previous computations has been made possible by the 
implementation of an automatic method for generating computer code for correlators of operators 
in various representations of the lattice rotational group starting from some basic operators. With 
this technique, inserting extra operators in our variational basis becomes an almost trivial task 
and the only limitations to the number of operators that can be used in the calculation are given 
by the available computational resources. In particular, using our method we are able to build 
a variational basis that includes torelon (i.e. states coupling to two Polyakov loops, which have 
an infinite mass in the infinite volume limit) and scattering operators; a variational ansatz on 
such a comprehensive basis allows us to disentangle genuine single-particle states from spurious or 
multi-particle resonances. The investigation reported in this work has been performed at a single 
value of the lattice spacing that has been fixed to a common value across the various SU(iV) gauge 
groups simulating at the gauge coupling corresponding to the deconfincment temperature when the 
temporal extension of the lattice is Nt = 6. The extrapolation to the continuum limit will be the 
subject of future investigations. 

The rest of the paper is organised as follows. The lattice setup will be briefly discussed in 
Sect. g. In Sect. || we shall review the variational method for extracting glueball masses, discuss 
scattering and torelon states, and briefly present our algorithm for automatically generating glueball 
operators in an irreducible representation of the rotational group. Our numerical results will be 
discussed in Sect. |[ with the extrapolation to large N given in Sect. ||. Finally, a summary of our 
findings together with an outlook on future developments will be the subject of Sect. g. 

2. SU(iV) gauge theories on the lattice 

The lattice discretisation of SU(iV) Yang-Mills theory used throughout this work is entirely conven- 
tional. We consider the system defined on an isotropic four-dimensional torus of linear size L. If a 
is the lattice spacing, the number of points in each direction is given by Nl = L/a. The volume of 
the system in physical units is given by V — L 4 . Directions are indicated with Greek symbols and 
lattice points are labelled with Latin letters. The fundamental degrees of freedom of the theory are 
the link variables U^i) £ SU(N). 11^(1) is associated to the bond (i; fi) that joins the site i with 
i + \x. The Wilson action for the lattice theory is given by 



S = P E (l-^ReTr(lVW)) , (2.1) 



t,fi,>h' 



where Ufj,v(i) is the parallel transport of the link variables along the elementary lattice plaquette 
stemming from point i and identified by the positive pair of directions (p, v). (3 is defined as 
j3 = 2iV/<7o, with go the bare gauge coupling. The Euclidean path integral reads 

(VU)e- s , (2.2) 

with VU the product of the Haar measures associated to each link. 

The observables we shall study in this work are glueball masses. For the continuum theory, 
at fixed quantum numbers J and for fixed excitation in the spectrum, the corresponding mass 
has corrections that can be expressed in a power series in \/N 2 . This is also true in the lattice 
strong coupling Q and has been shown to apply also for intermediate couplings [|8[ [46|, [l7|. In 
this work, we study the large N limit for fixed lattice spacing a as N varies. The lattice spacing is 
chosen in the domain where the behaviour of the discretised theory is dominated by the physics in 
the continuum limit, i.e. in a regime in which lattice corrections to continuum values of observables 
are modest. 

Since like QCD SU(iV) gauge theories arc defined by one dynamically generated scale, the 
physical value of the lattice spacing can be determined by measuring this length in units of a and 
assigning to the dynamical scale its continuum value. Conventionally, the square root of the string 
tension a is chosen to set the scale; fixing another quantity as a function of a amounts to corrections 
of order a 2 a to physical observables. In order to compare quantities at fixed lattice spacing across 
different SU(7V) groups, it proves useful to set the scale using the (pseudo)-critical coupling at 
fixed temporal extent. The justification for this choice (already successfully used in Q, |49| ) is that 
this quantity can be determined to a very high degree of accuracy. We adopt the deconfinement 
temperature at Nt — 6 to fix the lattice spacing to the same value when changing N. The choice 
of using /3 C at this particular value of Nj~ to set the scale is a compromise between the requirement 
of being in the continuum scaling regime and the practical convenience of having sufficiently large 
lattices in physical units for a computationally bearable value of the number of lattice sites N£ . In 
fact, Nl = 12 for (3 = (3 c (Nt — 6) gives a glueball spectrum in the scaling region and free from 
large finite size artefacts | J38[ . 

3. Extracting glueball masses 

In this section we present the construction of our operators and we review the general methodology 
for extracting glueball masses. While the standard variational procedure and the construction of 
operators in irreducible representations of the cubic lattice group is well known, this is, up to 
our knowledge, the first systematic attempt of inserting scattering and torelon operators into the 
variational set, in order to rule out from the spectrum contributions of these spurious states. 

Symmetries of the lattice spectrum 

At finite lattice spacing a the continuum rotation group is not an exact symmetry of the system. 
The full continuum rotational symmetry is dynamically restored only when a — > 0. On the lattice, 
eigenstates of the Hamiltonian have to fall into the irreducible representations of the octahedral 
point group Qo, the symmetry group of the cube. The octahedral point group has 5 irreducible 
representations A\, Ai, E, T\ and Ti respectively with dimensions 1, 1, 2, 3, 3. 
Since we are interested in the glueball spectrum of the gauge theory in the continuum, we need 
to consider Qo as a subgroup of the complete rotation group SO(3): irreducible representations 
of SO (3) are decomposed in terms of those of Qo- Irreducible representations of integer spin J in 
SO(3) restricted to Qo are referred to as subduced representations J J, Qo- When considered as a 
representation of Qo, the (2 J + 1) degeneracy of the continuum spin J state is split onto different 
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Figure 1: Set of basic prototypical paths used to construct operators in all the 20 R symmetry channels. 



irreducible representations of Qo ■ A simple example of this kind of pattern is the spin 2 (tensor) 
glueball, whose 5 polarisations are seen on the lattice as different states, 2 in the E and 3 in the 
T 2 representation of Qo . Due to the breaking of continuum rotational symmetry on the lattice, the 
aforementioned pattern of degeneracies is exact in the limit a — > 0, but it is only approximate at 
finite a. Comparing the measured glueball spectrum with the expected pattern of degeneracy can 
give information on the relevance of lattice artifacts. 

Near the continuum limit, it is possible to identify the masses of spin J glueballs by matching the 
patterns of degeneracies of the subduced representations J I Go from the degeneracy coefficients. 
We report these coefficients up to J = 4 in Tab. pi. 

For any given operator O on the lattice, we define a rota- 
tion transformation as TZi(O) where the index i labels all the 
elements of the group Qo- 

Since a generic representation of the group will not be irre- 
ducible, in order to create states that transform only in a given 
symmetry channel, we will need to create an appropriate linear 
combination of the rotations of the original operator. We define 
an operator in the irreducible representation R as 
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Table 1: Subduced representations 
J 4 Qo of the octahedral group up 
to J = 4. This table illustrates the 
spin content of the irreducible rep- 
resentations of Qo in terms of the 
continuum J. 



*<*>(t) = 5>f^(<5(t)) 



(3.1) 



„(«) 



The coefficients c\ appearing in the summation are obtained 
from the unitary operator that implements the change of basis 
from our choice of the representation in 24 dimensions into an orthonormal basis for each of the 5 
invariant subspaces [ p0| . 

By adding parity and charge conjugation to the group of pure rotations, we get the full symmetry 
group of glueball states on the lattice, which is referred to as Qq C - The group Qq C has a total 
of 20 irreducible representations labelled by R PC , where R indicates one of the 5 irreducible rep- 
resentations of Qo, P is the parity eigenvalue and C is the charge conjugation eigenvalue. As we 
will explain in details below, operators that are eigenstates of irreducible representations of Qo are 
constructed as traces of products of link variables along closed paths. Hence, eigenstates of C are 
given by the real (C = +1) or the imaginary (C = —1) part of the operator and adding and sub- 
tracting parity transformed operators gives definite -P states (respectively P = +1 and P = — 1). 
As usual we will label as R PC the ground state in a given symmetry channel, and its excitations as 
R PC followed by one or more *. 



Construction of the operators 

At finite volume, the single-particle glueball spectrum receives non-negligible corrections from 
multi-glueball states. Moreover, when the system is closed with periodic boundary conditions 
(like in our case) topological excitations wrapping the compact direction (torelons) with the same 
quantum numbers of glueballs appear; if not correctly accounted for, these states can affect sig- 
nificantly the measured glueball spectrum. In order to control these spurious contributions, we 
included in the variational set operators that best overlap with two-glueballs and torelon states. 
The variational procedure described in the next section will allow us not only to extract some 
excited states together with the lowest-lying state in each channel, but also to disentangle the 
contributions of spurious states. 

Before discussing the details of the variational procedure, we will describe the kind of operators 
used to define the variational set. 

All the operators described in the following are to be understood as gauge-invariant, vacuum- 
subtracted operators: 

6{t) = O(t)-(0\O(t)\0) . (3.2) 

Glueball operators 

The single-trace operator that we use to project onto glue- 
ball states is simply defined as 

0{t) = <j>{t) , (3.3) 

where <f>(i) is a zero-momentum operator given by the wall 
average over the temporal slice 
Table 2: Number of different glueball 

operators calculated in each of the 20 <£(*) = J- Y" #E, i) 4>(x,t) = Tx Y\ &i*(*) • 

symmetry channels. N L xeAs (l;A)eC 

(3.4) 
In our definition of the variational set we used a wide range 
of different closed loops C, with lengths ranging from 4 to 8 lattice spacings. The basic shapes we 
used in our calculations are summarised in Fig. |l|. 
In a given symmetry channel, we then correlate glueball operators of the form 

(3.5) 
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where the last term in second line is different from zero only when R = Aj . The number of 
different operators that we included in our variational set is summarised in Tab. for all the 
symmetry channels. 

Scattering operators 

An operator that projects onto scattering states of two glueballs is a double-trace operator. A trial 
operator that overlaps mainly onto scattering states has the form 

o(t) = m) - m))f , (3.6) 




Figure 2: Paths used for the construction of operators coupling with torelon states. Periodic Boundary 
Conditions apply at the edges represented by the dashed lines. 



where </>(£) is a zero-momentum operator defined as in Eq. ( |3.4| ). The operator appearing in the 
correlator matrix for the scattering states can be written using Eq. (|3.1|) as 



(3.7) 
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As in the previous case the last term vanishes for R ^ A^" 1 ", but the local subtraction of (<j>(t)) 

in the first term will appear in all the representations and is crucial in order to obtain the correct 

two-point function, even though (<f>(t)) alone, once appropriately symmetrised, would be different 

from zero only in the A{ channel. 

The local operators <f>(x, t) have been chosen among the closed loops in Fig. p] to create a variational 

set for the scattering states as big as the one used for glueball states. 

A technical aside: the aforementioned local subtraction in the scattering operators leads to the 

necessity of having access to the values of their vacuum expectation value (VEV) during the 

simulation, which implies that the VEVs must be calculated prior to the evaluation of the scattering 

operators. Moreover the VEV values have to be known with an accuracy higher than the most 

precise scattering operator. 

Torelon operators 

A torelon state is an excitation winding around the 
toroidal lattice in the spatial direction. The mass of a 
torelon state scales linearly with the lattice size L and 
the contribution of these states in the spectrum can be 
easily identified with a finite volume study. However, we 
can include in our variational set operators that mainly 
project onto torelon states and identify them on a single 
volume. 

Torelon excitations transform non-trivially under the cen- 
tre of the gauge group, and are characterised by their 
charge under this transformation (n-ality) . Since the glue- 
balls transform trivially under the centre of the gauge group, they can only couple to states that 
have zero n-ality. For this reason the torelon operators have been created from products of two 
Polyakov loops l v winding around opposite directions. We started from defining 
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Table 3: Number of different torelon op- 
erators calculated in each of the 20 sym- 
metry channels. 



°(*) = ^EE^^i^^) 



(3.8) 
where the sum over /j, runs on the spatial directions orthogonal to the one of the loops, and then 



the operator & R \t) can be constructed as in the previous cases by using Eq. (3.1). 
By choosing different shapes for the combination l v (x,t)l],{x + fia,t), we can obtain a fairly large 
variational set in all the 20 symmetry channels (see Tab. g). The shapes we used are pictorially 
shown in Fig. [2[ 



Blocking and smearing 

Masses of a lattice state R PC can be extracted from the behaviour of the correlator 



C (R \t) = (0|$ (fl)t (i)$ (fl) (0)|0> 



2 



AR)+, (3-9) 



= J2 <G|$ (R) |0> expC-m^t) 

G 

where G labels a generic eigenstatc of the action in an appropriate orthonormal basis. At large 
temporal distance, when only the lowest-lying state, labelled by the subscript 0, is contributing 
to the sum, it becomes possible to determine the mass of this state from the asymptotic single 
exponential decay 

lim C {R) (t) = |c | 2 exp(-m^ ) i) , (3.10) 

t—too 

where we have defined the overlap of the operator $^ on the state G as 



M 2 



(G|$ w |0) 2 (3.11) 



and the normalisation is chosen in such a way that 

Em 2 =1 • ( 3 - 12 ) 



In principle, Eq. (3.10Q allows one to extract glueball masses from the large t behaviour of 



correlators. However, it is already well known fed, BlL M, B2| that a good signal-over-noise ratio for 
the decay of this gluonic correlator can be obtained only at small temporal distances. Therefore it is 
mandatory to obtain the asymptotic single-exponential behaviour for t as small as possible. In fact, 
the single-exponential behaviour could be obtained at any t if there were only a single propagating 
state. The main idea underlying the variational technique is to enhance the propagation of a single 



2 



given state G in the sum of Eq. (3.9) and be able to find a normalised overlap \cq\ of order 1. To 



achieve such a result, standard smearing and blocking techniques [50, plfl have been implemented 
in the past, together with variational procedures. 

In this work, we used the improved blocking algorithm of [ p8[ . The basic idea is to create a 
configuration of super-links, obtained from an appropriate composition of blocking and smearing, 
from which we then build our "blocked" operators. The operator resulting from this procedure has 
a typical "size" given by the level of iterations of the algorithm. Blocked operators have a smaller 
sensitivity to fluctuations at the UV scale, and their overlap with the high-energy modes of the 
spectrum is reduced. 

In using the algorithm of |3q], we fixed the parameters of the algorithm to (p a ,Pd) — (0.40, 0.16) for 
the smearing part (see Fig. pi) , while a blocked link is obtained multiplying two consecutive smeared 
links. This choice has been shown to produce operators with |cg| ~ 0(1) for all symmetry chan- 
nels. 

As a result, operators constructed at different blocking levels Nf, will allow us to enlarge the vari- 
ational set introduced above. For example, each of the glueball operators of Tab. |] can be built 
using blocked super-links at different levels and, in our variational procedure, we were thus able to 
use a set that is effectively Nj, = 4 times bigger than the original one. 

A remark is in order here. The blocking procedure may lead to the construction of gauge variant 
operators if the size of the lattice along the blocked direction is not multiple of 2 Nb ~ 1 and one 
is measuring operators that wind along this direction. This problem can be circumvented mixing 
different lower blocking levels [p4[. 
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Figure 3: Example of smearing of a single link. Sum of staples and diagonal staples are weighted with 
two independent parameters, respectively p a and pd- 



Variational procedure 

Let us briefly describe the procedure we have used to obtain the mass of a lattice state starting from 
a variational basis constructed as discussed above. In the following, we denote by $„ an operator 
that creates a state in the R channel, just as defined in Eq. (3.1), and the index a £ 1, . . . , N labels 



one of the N operators in the variational set. Since we will be dealing with one channel at a time, 

from now on, we will drop the (R) superscript for the sake of clarity. 

The variational ansatz consists in finding an appropriate linear combination of the basis operators 



*( f ) = S ««*«(*) 



(3.13) 



whose correlator contains a single propagating state for which the mass is minimised by the optimal 
v a . The optimisation of v a can be solved as a generalised eigenvalue problem. We start by measuring 
a correlator matrix between all the possible choices of operators in the variational set 



C ap {t) = ^(0|$i(< + r)$^(r) 



(3.14) 



Then we diagonalise the matrix C(t) = C 1 (0)C(t) to solve the eigenvalue equation 



(3.15) 



where we choose t = 1. Each eigenvalue gives the mass of a state in the chosen symmetry channel, 
thus the largest eigenvalue corresponds to the ground state and its eigenvector is used to construct 



the operator that best overlaps onto it following Eq. ( 3.13 ). In the present work we studied the 
eigenvectors corresponding to the 5 highest eigenvalues to construct operators $(£)», i — 1, • • • , 5, 
whose diagonal correlators Cu (£) are fitted with the single-cosh ansatz 



Cu{t) = \c l \ 2 cosh.{m l t~N L /2) 



(3.16) 



which assumes that the exponential of Eq (3.G) is dominated by only one state with mass m.;. 
This procedure allows us to extract also masses of excited states since the operators $,(£) that 
decay with a mass TOj ^ m,Q within errors are generally good approximations of excited glueballs 
wavefunctions. 

A large variational set can result in a singular correlation matrix if accidental linear dependencies 
between the operators develop. To prevent singularities in the correlation matrix, we have denned a 
posteriori the variational basis in each channel as the largest set of operators for which the inversion 
of C can be performed with a numerical accuracy of order 10 -13 or better. 

From the point of view of the implementation, it should be clear that the simulation code must 
be able to handle a huge amount of different shapes for the construction of the operators. Moreover, 
it is not known a priori which are the relevant shapes in the definition of a good variational set. 
Thus the simulation should be flexible enough to allow us to change the used shapes with little effort. 
In our project we have decided to develop a totally automated environment that is able, at compile 
time, to create the code to measure all the possible shapes in all the possible symmetry channels, 



with the only input of a string representing each shape we want to include in the calculation. This 
environment has been created by means of a Mathematica code that, starting from a set of shapes 
and symmetry channels, writes the plain functions to measure the properly symmetrized operators 
and their correlators. Such a solution turns out to be very efficient, since all the running time is 
used only in measuring the operators. In the end our environment allows us to use a large basis of 
operators, whose size is limited only by the amount of available memory and computer time, and 
that can hence be increased at will with no programming effort. 

4. Numerical results 

Our calculations are performed for 3 < N < 8. As discussed in Sect, g, at each N we simulate at 
the coupling for which the system is at the deconfinement temperature on a Nt — 6 lattice. For 
N = 3,4,6,8, (3 c (Nt = 6) is available from previous calculations |J2l |5^, p3. For the remaining 
values of N, we have computed (3 c (Nt) as follows p2|, |53"[ . On a lattice with geometry jV£ x Nt we 
computed first /3 C for Nt — 5 using an extrapolation to the thermodynamic limit of the position of 
the peak of the Polyakov loop susceptibility and then we assumed the coefficient of the l/iV£ term 
from that fit to extrapolate to the thermodynamic limit the peak position of the Polyakov loop 
susceptibility from results obtained at Nt = 6 at one value of Nt in the asymptotic regime. Our 
analysis confirms the expectation that for N — 5 and N = 7 the deconfinement phase transition is 
first order; the details of this calculation will be reported elsewhere. Tab. summarises the relevant 
values of (3 c (Nt — 6) at the values of iV" used in our investigation. The values of (3 c (Nt = 6) at 
N = 3,4,6,8 are taken from Ref. |53] (Ref. [E34] quotes compatible results), while the values at 
N = 5,7 have been computed as a part of this work. 

The parameters used in the simulations are also summarised in Tab. m A typical run has N T 
thermalisation sweeps after which we start measuring the correlators. A Monte Carlo step consists 
of a compound sweep in which one heat -bath update is followed by 4 over-relaxation sweeps. After 
the thermalisation process we perform Nmc sweeps. We chose to measure every N compoun d sweeps 
to reduce autocorrelation between the measures. The N measures measure sweeps are further divided 
in Nbi ns bins; each bin is an average over N W idth measures. The total set of measures to be analysed 
is then Nuns x runs, where each run is independent of all the others. 
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17.1068(30) 
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100k 
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Table 4: Values of the critical couplings for the deconfinement temperature at Nt — 6 and parameters of 
the Monte Carlo simulations on lattices with N L points for 3 < N < 8. 

After performing the variational calculation, the diagonal elements of the correlation matrix 



are fitted with the single-cosh ansatz of Eq. ( 3.16 ), which assumes that only one state dominates 



the signal. Unitarity imposes \ci\ 2 < 1. In practise, we are often able to obtain overlaps of the 
order of 0.95, which proves the validity of the original variational ansatz. As a consequence, the fit 
generally works very well on the range 1 < t < 4. 

In our investigation, we start by excluding at first the scattering states from our variational 
basis. Our numerical results are summarised in Tabs. |5jT0, where we list for all the gauge groups 



the channels for which we were able to extract at least the groundstate mass; where more than a 



mass has been extracted, excitations are denoted with a number of * indicating their position in 
the spectrum, according to the standard convention. For each gauge group, we are able to obtain a 
signal in most of the J PC channels, with a good quality of the fit. We have highlighted in boldface 
the states for which the x 2 is bigger than 1.25; similarly, we have indicated in boldface states for 
which the overlap is less than 0.85 (possibly indicating that the single-cosh fit might be a sub- 
optimal approximation) or statistically bigger than 1 (the error on the overlap coming from the 
fit, which is not reported, is at most of the order of 0.1). Only the states that are free from those 
problems (which are a large majority) will be used for the large N extrapolation. 

The symmetry channel for which we are able to extract more states is the Aj . Symmetry 
channels whose groundstate is the same state in the continuum (e.g. the E ++ and the T 2 ++ , whose 
lowest-lying state is expected to be the 2 ++ glueball) have compatible groundstate mass with the 
expected degeneracy at our lattice spacing, confirming that for our choice of parameters we are in 
the scaling region. 

The states obtained after the variational procedure can be decomposed into their projection 
onto the pure glueball states and onto the torelons: 

$ = a G f G + otT$T , (4.1) 

where $g and $r are normalised to unity. This allows us to define the pure glueball component 
and the torelon component of those states respectively as 

IqgI 2 Iq^tI 2 

miXG = i — 12T1 — [2 and miXT = i — im — [2 • ( 4 - 2 ) 

|ckg| + larr \oig\ +\aTr 

For all gauge groups, there is a high mixing between narrow glueball trial states and torelon states 
in the first excitation of the E ++ and in the second excitation of the A^ + . Other states with a 
consistent mixing with the torelons are the 1% and the Tj , the latter mostly for N = 3,4. 
In order to guide the eye, in Tabs. |5|-|l0| we have labelled by (2T) the states that have an overlap 
with the torelons of order 0.35 or bigger and by (2T?) the states with an overlap between 0.15 and 
0.35. For all other states, the overlap with the torelons is expected not to influence significantly 
the large- volume scaling of their mass. 

We now turn to the problem of mixing of narrow resonances with scattering states. Since a 
calculation involving scattering states is much more demanding in terms of computer time, we use 
the results from the computation involving only single-particle and torelon operators to target the 
channels where mixing with mult i- particle states is expected to affect significantly the results. At 
large N, this is expected to happen for the excited states that are close to twice the energy of the 
groundstate. It is then clear that the channel in which scattering states can potentially influence 
the measured spectrum in a relevant way is the A± , where we can extract several excitations. 

In our investigation, we limit the study of scattering states to the Ai channel, for which we 
perform a calculation on a larger set of operators, including multi-glueball operators in addition to 
the single- glueball and torelon operators used before. In fact, we perform calculations on separate 
sets of operators (the full set and the sets obtained excluding in turn scattering, torelon and single- 
glueball operators). Our results are reported in Tabs. [ll||Lq . In analogy with the calculation 
involving only single-particle and torelon operators, we define the mixing coefficients mixc, mixg 
and mixT that give the projection of the variational states respectively onto the single-glueball, 
scattering and torelon subspaces. The mass spectrum obtained with the different choices for the 
operator sets is illustrated in Figs. &H. The remarkable property shown by this calculation is 
that when only scattering and torelon operators are used the lowest-lying state has a mass that is 
much higher (roughly by a factor of two) than the mass of the groundstate extracted with the full 
variational basis. Moreover, the latter appears always when single-particle operators are included 
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in the calculation. This is an indication that our multi-glueball set of operators project only on 
scattering states, as it should be. The scattering state seems to be slightly above the first excited 
single-glueball excitation at any value of N, and corresponds to the state identified as Af + ** in 
the calculation that excludes the scattering operators (note the relatively large overlap of this state 
with the torelon operators in that calculation). Another feature we notice is that the number of 
excitations we are able to extract depends on the subset of operators used. This is hardly surprising, 
given the variational nature of the procedure. A similar analysis in the A^ + channel shows that 
the A^ state identified with a variational calculation excluding scattering operators is in fact a 
scattering state. 

We performed an additional check of finite volume effects in the A^ + channel studying SU(3) at 
Nl = 18, with the results illustrated in Fig. H. This test shows the expected result that torelons are 
not visible anymore in the lowest-lying spectrum on this larger volume. In addition, the computed 
mass spectrum in the -AJ" channel is compatible with the one measured for Nl = 12. Since finite 
size effects are expected to be more relevant for the lightest glueballs, this check indicates that finite 
size corrections are negligible for Nl = 12. 

In conclusion, our variational set enables us to correctly disentangle single-particle states from 
scattering states and torelon excitations. Moreover, the mixing coefficients give a measure of the 
contamination of the single-particle states from unwanted contributions. In the following, we con- 
sider single-glueball states all the states with a projection onto the single-glueball operator sector 
of 85% or larger and we disregard states with a projection of 20% or larger onto the torelon sector. 
Scattering states are identified as the states in the spectrum with a significant (30% or larger) 
projection onto the scattering sector that are degenerate in mass with the lowest-lying state in the 
S + T basis when the latter has the largest component in the scattering sector. 

5. The large N spectrum 



In [Bq, 38 it has been shown that the lowest-lying states of the spectrum can be extrapolated to 



the large N limit using only the leading large N correction: 

am G (N) = am G (oo) + c/N 2 , (5.1) 

where c is a coefficient of order one which depends on the symmetry channel. We perform a similar 
extrapolation for the states extracted in our work. We include in the fit all the states we have 
extracted from the variational procedure, indicating in boldface in the tables those who have been 
identified as potentially problematic in the previous section because of the overlap or of the larger 
X 2 . The obtained large N spectrum is reported in Tabs. [lTJ-glp. Generally, when a sufficient 
number of reliable states is present, the inclusion of the more dubious states does not change the 
result within errors. This is shown explicitly in Tab. E3, where the extrapolations with and without 
the less reliable states are compared. Both extrapolations are displayed in Figs. [10-19 for all the 



quantum numbers for which a large N limit can be extracted. We find that the ansatz (5.1) works 
for all the measured states (including the excitations) for N > 3. In general, the central value of c 
is found to be small (always of order one or below), as it is expected for a generic coefficient in a 
well-behaved expansion. For most of the states we find only modest corrections to the N = oo value 
of the mass: with a few exceptions, c is compatible with zero and a fit with only the leading term 
amcioo) in Eq. (|5.l|) gives a result that is compatible with the fit that includes also the 0(1/N 2 ) 
correction. In the following, we shall use the fit results obtained with the ansatz (|5.l|). 



1 Note that for the A~[ channel we have indicated with A^ + (S) the scattering states identified in the previous 
section and relabelled genuine single-particle excitations excluding these states from the spectrum. An analogous 
relabelling has been done in the A^~ channel. 
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As expected, the channel for which we extract the highest number of states is the Af + , where 
we clearly see the groundstate, the first two excitations and a scattering state. The large N 
extrapolation is shown in Fig. |10L The groundstate is clearly visible also in the A± channel, lying 
slightly below the first excitation of the Af + channel (Fig. 111]). Fig. O] shows the extrapolation of 
the groundstate in the A 2 channel, which is the only state we can extract for the A<i irreducible 
representation of the octahedral group. It is also clearly shown the effects of including less reliable 
masses in the extrapolation fit. In the E channels, we get the groundstate and the first excitation 
for the E ++ (Fig. [R|), in addition to a torelon state that we discard because it is not a spectral 
state. The ground state and the first excitation are also visible in the E h (Fig. [u]), while only the 
groundstate is present in our extrapolated data for the E^ and E~~ (Fig. |l5|). The groundstate 
is clearly visible in the T^ + and T 1 ++ channels (Fig. Hq), while also the first excitation can be 
seen in the T^ channel (Fig. |l7| ). Finally, two states appear in the T 2 + channel (Fig. |l^), with 
the groundstate determined for the T 2 + , ^2 anc ^ the T 2 channels (Fig. |l^). The large N 
extrapolation preserves the degeneracy of states that have the same spin in the continuum limit. 
The only exception to this general fact is observed for the (E + ~ ,T 2 ) pair. We notice however 
that the extrapolation to N — > 00 of the T 2 groundstate only relies on two values at finite N. 
This probably explains the splitting between the two states at N — 00. 

We have also extrapolated to the large N limit the scattering states identified in the A^ + and 
A^ + channels. For the former state, for which we have better control, the extrapolated value of the 
energy is compatible with twice the mass of the groundstate, in agreement with the expectation 
that particles do not interact at large N. The energy of the Aj" + (S) is slightly smaller than twice 
the mass of the corresponding groundstate, but this should probably be ascribed to the lack of 
control over the extrapolation to the large N limit, which had to rely only on two finite- N values. 

The single-glueball spectrum determined in this work is plotted in Fig. Gu and it is compared 
with the known spectrum at the same lattice spacing taken from Ref. fl3q|. While the latter work 
achieve a better precision for the A^ + , the E ++ and the A^ + * , in this study we are able to 
measure seventeen more states. Moreover, our results for the common states are compatible with 
those reported in fl38|| . 

The determination of masses in various symmetry channels opens the possibility to study Regge 
trajectories at large N. A mandatory first step in this context is the identification of the continuum 
spin corresponding to each lattice state. The identification of the continuum spin is in general a 
very complicated task (see e.g. Ref. p2| for a discussion). In this work we use the simple-minded 
approach of looking at the subduced representations reported in Tab. [fl, with the appearance of 
the same continuum state in various Qo representations taken into account whenever the expected 
degeneracies are unambiguous. We associate all states in the A\ representation with spin J = 0, the 
T\ states are associated with J = 1 and the only J = 3 state is obtained from the A-i lattice state. 
For J = 2, we average the E and the corresponding T2 states when possible, except for J PC = 2 H 
for which we take only the E^ , having a poor control over the errors from the T 2 extrapolation. 
The resulting Chew-Frautschi plot is shown in Fig. Ell In order to obtain the string tension in lattice 
units at our lattice spacing, we have performed a large N extrapolation of the data for T c / \fa in 
Refs. H, 0' from which we obtain T c /^ = 0.5851(32) for N T = 6 and N = 00. Although our 
results are in qualitative agreement with Ref. J42|, in order to check reliably predictions coming 
from various models more excitations and mostly higher spin states would be needed. 

6. Conclusions 

In this work, we have studied numerically on the lattice the glueball spectrum in Yang-Mills SU(iV) 
gauge theories in the large N limit. Using an automated technique for constructing trial wave 
functionals in all possible symmetry channels, we have built a large variational basis that has 
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enabled us to obtain a large number of states, including some excitations. Moreover, the inclusion 
of functionals that best overlap with scattering and torelon states has allowed us to unambiguously 
exclude multi-particle states or finite-size artefacts from the spectrum of narrow resonances. This 
is a significant advance in our understanding of the large N glueball spectrum from first principles. 
Our calculation being at fixed lattice spacing, a natural development of this work would be the 
determination of the spectrum in the continuum limit. We expect that as the lattice spacing 
becomes finer, more and more states will be pinned down by our variational calculation. In this 
respect, it might be helpful to use anisotropic lattices, like in Ref. f32|| . Together with a correct 
identification of the continuum spin states corresponding to the glueball states classified according to 
the irreducible representations of Go > a similar calculation would open the way to a comprehensive 
investigation of Regge trajectories, which are expected to be a more clean signature of the spectrum 
at large N than at N — 3. With little or no modification, the technique we have presented in this 
work will also prove helpful in related problems, like the lattice study of glueballs in QCD and the 
study of the low-energy spectrum of confining flux tubes. 
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5.8945 and N L = 12 
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1.324(42) 


0.98 


0.61 


0.9828 


0.0172 


J 2 


1.79(11) 


0.97 


0.16 


0.9635 


0.0365 


J 2 


1.79(11) 


0.96 


0.38 


0.9887 


0.0113 


T+-(2T) 


1.87(14) 


0.91 


0.33 


0.5195 


0.4805 


T 2 +_ (2T?) 


1.85(13) 


0.89 


1.00 


0.7641 


0.2359 


T 2 +~(2T?) 


2.13(19) 


1.12 


0.19 


0.6719 


0.3281 


r 2 ~ 


2.26(28) 


0.97 


0.18 


0.9968 


0.0032 


J 2 


1.665(93) 


0.94 


0.73 


0.9803 


0.0197 


J 2 


1.614(82) 


0.89 


0.01 


0.9788 


0.0212 


J 2 


2.13(26) 


0.88 


0.70 


0.9722 


0.0278 



Table 5: The measured spectrum of pure SU(3) gauge theory from lattice simulations at the parameters 
shown in the header. 
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SU(4) at /3 = 


10.789 and N L 


= 12 




R PC 


amo 


M 2 


X 2 


mixc 


mixT 


Af + 


0.821(15) 


0.99 


1.03 


0.9928 


0.0072 


A++* 


1.415(47) 


0.96 


1.68 


0.9871 


0.0129 


A++*(2T?) 


1.596(73) 


0.97 


0.74 


0.7871 


0.2129 


A+ + ** 


1.97(18) 


0.85 


0.43 


0.9842 


0.0158 


Af- 


2.79(86) 


0.86 


0.04 


1.0000 


0.0000 


AI + 


1.423(54) 


0.89 


0.21 


0.9991 


0.0009 


A^ + * 


2.67(63) 


1.08 


0.67 


0.9927 


0.0073 


A++(2T) 


1.83(15) 


0.85 


1.55 


0.0271 


0.9729 


A++*(2T?) 


2.66(66) 


0.87 


0.11 


0.8390 


0.1610 


At~ 


2.08(19) 


0.99 


0.27 


0.9998 


0.0002 


Ar 


2.34(37) 


0.87 


0.47 


0.9291 


0.0709 


E++ 


1.270(41) 


0.94 


0.01 


0.9787 


0.0213 


E++ 


1.346(47) 


1.00 


0.46 


0.9553 


0.0447 


E++{2T) 


1.634(75) 


1.03 


0.43 


0.3336 


0.6664 


E++{2T) 


1.711(91) 


1.07 


0.01 


0.2347 


0.7653 


E+- 


2.85(70) 


1.17 


0.66 


0.9728 


0.0272 


E~+ 


1.732(97) 


0.95 


0.03 


0.9922 


0.0078 


E~+ 


1.653(92) 


0.84 


0.48 


0.9953 


0.0047 


T++ 


2.13(22) 


0.88 


0.35 


0.9911 


0.0089 


Tf- 


1.638(83) 


0.99 


1.19 


0.9867 


0.0133 


T+- 


1.639(77) 


0.98 


0.20 


0.9634 


0.0366 


T+- 


1.718(88) 


1.03 


0.01 


0.9657 


0.0343 


T+" 


1.79(13) 


0.85 


0.05 


0.9464 


0.0536 


T+-* 


2.15(19) 


1.15 


0.02 


0.9522 


0.0478 


Tf - (2T) 


2.17(28) 


0.88 


0.45 


0.5640 


0.4360 


Tf - (2T) 


2.45(39) 


0.99 


0.02 


0.4546 


0.5454 


IT+ 


2.37(36) 


0.90 


0.02 


0.9842 


0.0158 


J 2 


1.333(47) 


0.99 


0.01 


0.9858 


0.0142 


1 2 


1.76(11) 


0.89 


0.17 


0.9910 


0.0090 


T+-(2T) 


1.96(16) 


0.90 


0.48 


0.6320 


0.3680 


T+-(2T) 


2.24(23) 


1.16 


0.20 


0.5183 


0.4817 


T+-(2T) 


2.12(20) 


1.00 


0.22 


0.2399 


0.7601 


r 2 ~ 


2.21(25) 


0.89 


1.75 


0.9883 


0.0117 


1 2 


1.609(87) 


0.86 


0.04 


0.9842 


0.0158 


1 2 


1.678(95) 


0.88 


0.01 


0.9597 


0.0403 



Table 6: The measured spectrum of pure SU(4) gauge theory from lattice simulations at the parameters 
shown in the header. 
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SU(5) at /3 = 


: 17.107 and N L = 12 




R PC 


arriQ 


M 2 


x 2 


mix<3 


mix<r 


At + 


0.798(14) 


0.98 


0.70 


0.9901 


0.0099 


A++* 


1.426(57) 


0.92 


0.01 


0.9864 


0.0136 


A++(2T?) 


1.675(99) 


0.94 


1.61 


0.7006 


0.2994 


A+ + ** 


2.03(20) 


0.88 


0.40 


0.9986 


0.0014 


AI + 


1.406(57) 


0.88 


0.02 


0.9975 


0.0025 


A7+* 


2.28(29) 


0.92 


0.48 


0.9877 


0.0123 


A++(2T) 


1.88(16) 


0.85 


0.08 


0.0750 


0.9250 


Ai- 


2.19(22) 


1.11 


0.35 


0.9980 


0.0020 


E++ 


1.300(42) 


0.95 


1.00 


0.9830 


0.0170 


E++{2T) 


1.628(87) 


0.94 


1.73 


0.2965 


0.7035 


E++{2T) 


1.597(84) 


0.89 


0.03 


0.5388 


0.4612 


E++* 


1.91(14) 


1.05 


0.20 


0.9619 


0.0381 


E— 


2.05(22) 


0.83 


0.01 


0.9737 


0.0263 


E~+ 


1.72(10) 


0.90 


0.01 


0.9905 


0.0095 


E-+* 


2.48(43) 


0.97 


0.05 


0.9956 


0.0044 


T++ 


2.09(21) 


0.91 


0.03 


0.9952 


0.0048 


T++ 


2.03(21) 


0.82 


0.30 


0.9961 


0.0039 


T+- 


1.598(80) 


0.95 


0.00 


0.9820 


0.0180 


T+- 


1.625(76) 


0.96 


0.82 


0.9570 


0.0430 


T+-* 


1.99(16) 


0.99 


0.11 


0.9107 


0.0893 


Tf - (2T?) 


2.16(28) 


0.81 


0.70 


0.8117 


0.1883 


ir + 


2.30(30) 


0.90 


1.19 


0.9953 


0.0047 


J 2 


1.261(41) 


0.90 


0.13 


0.9738 


0.0262 


1 2 


1.321(45) 


0.94 


0.03 


0.9784 


0.0216 


rp+ + 

1 2 


1.391(48) 


1.02 


0.19 


0.9931 


0.0069 


rp+ + * 

1 2 


1.76(12) 


0.86 


0.39 


0.9809 


0.0191 


rp+ + * 

1 2 


1.91(14) 


0.99 


0.77 


0.9798 


0.0202 


1 2 


1.90(15) 


0.84 


0.55 


0.9670 


0.0330 


T 2 +~(2T?) 


1.94(15) 


0.94 


0.03 


0.8434 


0.1566 


T+-(2T) 


2.03(20) 


0.87 


0.18 


0.2209 


0.7791 


T 2 - 


2.01(21) 


0.80 


0.65 


0.9955 


0.0045 


T -+ 

1 2 


1.71(10) 


0.92 


0.71 


0.9865 


0.0135 


rp-+* 

1 2 


1.93(12) 


1.13 


0.56 


0.9895 


0.0105 


rp h* 

1 2 


2.31(31) 


0.95 


0.08 


0.9601 


0.0399 


rp K 

1 2 


2.29(29) 


0.91 


0.86 


0.9490 


0.0510 



Table 7: The measured spectrum of pure SU(5) gauge theory from lattice simulations at the parameters 
shown in the header. 
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SU(6) at /3 = 


24.845 and N L 


= 12 




R PC 


airiQ 


|co| 2 


x 2 


mixG 


mixT 


At + 


0.785(14) 


0.96 


0.20 


0.9959 


0.0041 


A++(2T) 


1.572(83) 


0.82 


0.54 


0.5763 


0.4237 


A++* 


2.37(28) 


1.19 


0.05 


0.9842 


0.0158 


A-, + 


1.417(59) 


0.89 


1.51 


0.9996 


0.0004 


a + 2 - 


1.85(14) 


0.80 


0.74 


0.9993 


0.0007 


At~ 


2.59(51) 


0.99 


0.03 


0.9991 


0.0009 


E++ 


1.257(43) 


0.89 


0.40 


0.9728 


0.0272 


E++{2T) 


1.607(86) 


0.92 


0.83 


0.2433 


0.7567 


E++*(2T?) 


1.97(14) 


1.09 


0.16 


0.8404 


0.1596 


E— 


2.33(28) 


1.07 


0.02 


0.9913 


0.0087 


E~+ 


1.76(11) 


0.93 


0.75 


0.9929 


0.0071 


E~+ 


1.71(11) 


0.85 


1.16 


0.9951 


0.0049 


E -+* 


2.56(45) 


1.07 


0.01 


0.9886 


0.0114 


T++ 


2.30(28) 


0.99 


0.01 


0.9924 


0.0076 


T++ 


2.29(33) 


0.92 


0.66 


0.9919 


0.0081 


Tf- 


1.501(66) 


0.86 


0.45 


0.9860 


0.0140 


T+- 


1.658(77) 


0.97 


0.66 


0.9539 


0.0461 


T+-* 


1.94(16) 


0.93 


0.11 


0.9443 


0.0557 


r r + 


2.43(43) 


0.93 


0.03 


0.9859 


0.0141 


1 2 


1.263(39) 


0.90 


0.74 


0.9881 


0.0119 


1 1 


1.367(46) 


0.99 


1.96 


0.9841 


0.0159 


1 2 


1.99(14) 


1.05 


0.56 


0.9704 


0.0296 


1 2 


2.02(16) 


1.07 


1.71 


0.9879 


0.0121 


1 2 


1.98(16) 


0.93 


0.14 


0.8618 


0.1382 


T +- 
1 2 


2.08(18) 


1.02 


0.78 


0.9465 


0.0535 


T 2 - 


2.31(30) 


0.99 


0.15 


0.9972 


0.0028 


1 2 


1.85(12) 


1.05 


0.18 


0.9416 


0.0584 


1 2 


1.74(11) 


0.91 


0.13 


0.9765 


0.0235 


rp- + 

1 2 


1.63(10) 


0.80 


0.88 


0.9790 


0.0210 



Table 8: The measured spectrum of pure SU(6) gauge theory from lattice simulations at the parameters 
shown in the header. 
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k 


5U(7) at /3 = 


= 33.995 and N 


l = 12 




R PC 


a?no 


|co| 2 X 2 


mixc 


mixr 


At + 


0.820(15) 


0.98 0.43 


0.9976 


0.0024 


A++* 


1.564(61) 


1.08 1.08 


0.9920 


0.0080 


A++{2T) 


1.744(89) 


0.99 1.29 


0.6089 


0.3911 


A+ + ** 


2.30(27) 


1.06 0.07 


0.9970 


0.0030 


A-, + 


1.452(58) 


0.92 0.07 


0.9994 


0.0006 


A,+* 


2.28(32) 


0.89 0.23 


0.9968 


0.0032 


A++(2T) 


1.84(13) 


0.78 1.12 


0.0416 


0.9584 


At + 


2.18(25) 


0.88 0.44 


0.9231 


0.0769 


A\- 


2.43(43) 


0.84 0.23 


0.9996 


0.0004 


E++ 


1.296(44) 


0.95 0.03 


0.9841 


0.0159 


E++ 


1.303(46) 


0.93 1.31 


0.9723 


0.0277 


E++{2T) 


1.649(89) 


0.92 1.44 


0.2811 


0.7189 


E++{2T) 


1.72(10) 


0.93 1.19 


0.3631 


0.6369 


E++* 


1.81(12) 


0.91 0.73 


0.8825 


0.1175 


E+- 


2.58(62) 


0.87 0.39 


0.9948 


0.0052 


E+- 


2.63(64) 


0.85 0.04 


0.9950 


0.0050 


E— 


2.42(30) 


1.11 1.50 


0.9904 


0.0096 


E—(2T) 


2.79(66) 


1.05 0.23 


0.5190 


0.4810 


E~+ 


1.78(11) 


0.91 0.24 


0.9974 


0.0026 


E-+* 


2.59(42) 


1.09 0.66 


0.9958 


0.0042 


T++ 


2.27(25) 


0.99 1.92 


0.9962 


0.0038 


T++ 


2.42(31) 


1.15 0.44 


0.9874 


0.0126 


T+- 


1.543(73) 


0.89 0.11 


0.9628 


0.0372 


T+- 


1.617(82) 


0.93 0.52 


0.9581 


0.0419 


T+- 


1.596(79) 


0.90 0.15 


0.9734 


0.0266 


T+~* 


1.89(15) 


0.85 0.01 


0.9895 


0.0105 


T— 


2.12(22) 


0.87 0.60 


0.9334 


0.0666 


Tr 


2.46(33) 


1.11 0.88 


0.9548 


0.0452 


Tr 


2.33(37) 


0.84 0.07 


0.9092 


0.0908 


ir + 


2.48(46) 


0.99 1.98 


0.9605 


0.0395 


1 2 


1.322(45) 


0.95 0.24 


0.9867 


0.0133 


rp+ + 

1 2 


1.367(50) 


0.97 0.99 


0.9879 


0.0121 


rp+ + 

1 2 


1.437(50) 


1.04 1.73 


0.9799 


0.0201 


rp+ + * 

1 2 


1.96(15) 


0.99 0.21 


0.9780 


0.0220 


1 2 


2.07(17) 


0.98 1.27 


0.9682 


0.0318 


T+-(2T) 


2.05(18) 


0.93 1.73 


0.6236 


0.3764 


T+"(2T) 


2.23(23) 


1.04 0.21 


0.2859 


0.7141 


T 2 - 


2.45(33) 


1.13 0.20 


0.9886 


0.0114 


T -+ 

1 2 


1.77(11) 


0.97 0.40 


0.9738 


0.0262 


rp h* 

1 2 


2.36(34) 


0.98 1.11 


0.9618 


0.0382 



Table 9: The measured spectrum of pure SU(7) gauge theory from lattice simulations at the parameters 
shown in the header. 
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SU(8) at /3 = 


44.496 and JVj 


. = 12 




R PC 


airiQ 


M 2 


X 2 


mixG 


mixr 


Af + 


0.785(16) 


0.96 


0.77 


0.9993 


0.0007 


A++* 


1.408(63) 


0.90 


0.52 


0.9966 


0.0034 


A-, + 


1.401(62) 


0.83 


1.48 


0.9995 


0.0005 


At~ 


2.21(27) 


1.08 


2.38 


0.9981 


0.0019 


E++ 


1.339(53) 


0.96 


0.38 


0.9937 


0.0063 


E++{2T) 


1.78(12) 


0.94 


0.34 


0.5444 


0.4556 


E++* 


1.80(14) 


0.88 


0.44 


0.9330 


0.0670 


E+- 


2.67(72) 


0.93 


0.02 


0.9815 


0.0185 


E— 


2.24(30) 


0.92 


1.59 


0.9876 


0.0124 


E— 


2.50(41) 


1.15 


0.27 


0.9806 


0.0194 


E~+ 


1.90(14) 


1.05 


0.55 


0.9976 


0.0024 


T++ 


2.25(33) 


0.87 


0.32 


0.9915 


0.0085 


T+- 


1.693(97) 


1.00 


0.07 


0.9742 


0.0258 


T+" 


1.682(97) 


0.98 


1.28 


0.9847 


0.0153 


1+- 


1.590(99) 


0.88 


0.09 


0.9354 


0.0646 


T+~* 


1.98(17) 


0.97 


0.60 


0.9748 


0.0252 


T+-* 


2.10(24) 


1.04 


0.52 


0.9803 


0.0197 


Tf - (2T?) 


2.16(27) 


0.88 


0.70 


0.7766 


0.2234 


Tf-(2T?) 


2.47(43) 


1.14 


1.15 


0.8164 


0.1836 


rr + 


2.45(43) 


0.95 


0.80 


0.9909 


0.0091 


ir + 


2.57(57) 


0.98 


0.53 


0.9931 


0.0069 


J 2 


1.323(51) 


0.93 


0.00 


0.9853 


0.0147 


J 2 


1.451(56) 


1.05 


1.23 


0.9687 


0.0313 


J 2 


1.393(57) 


0.97 


0.39 


0.9749 


0.0251 


rrn+ + * 

J 2 


1.90(15) 


0.93 


1.67 


0.9878 


0.0122 


1 2 


2.15(21) 


1.12 


0.72 


0.9232 


0.0768 


1 2 


2.07(21) 


0.97 


0.02 


0.9754 


0.0246 


T +- 


1.99(18) 


0.89 


1.25 


0.9644 


0.0356 


T+-(2T) 


2.11(24) 


0.92 


1.41 


0.3123 


0.6877 


T+-(2T) 


2.26(28) 


1.01 


1.25 


0.3711 


0.6289 


r 2 — 


2.20(32) 


0.93 


0.54 


0.9870 


0.0130 


TF" 


2.12(31) 


0.82 


1.05 


0.9904 


0.0096 


1 2 


2.01(16) 


1.19 


0.02 


0.9786 


0.0214 


1 2 


1.79(14) 


0.95 


0.95 


0.9628 


0.0372 


1 2 


2.22(33) 


0.87 


0.41 


0.9306 


0.0694 



Table 10: The measured spectrum of pure SU(8) gauge theory from lattice simulations at the parameters 
shown in the header. 
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SU(3) at /3 = 5.8945 and N L = 


12 




R PC 


aTO(cr) 


\Cn\ 2 


x 2 


mix G 


mixg 


mix-r 


A++ 


0.792(14) 


0.98 


0.02 


0.9664 


0.0180 


0.0156 


A++* 


1.370(50) 


1.01 


0.01 


0.8242 


0.0870 


0.0888 


<++** 


1.609(97) 


0.90 


0.74 


0.8060 


0.1753 


0.0187 


A+ + ±±* 


1.86(13) 


0.99 


0.05 


0.9197 


0.0531 


0.0273 


At + 


0.792(14) 


0.98 


0.02 


0.9830 


0.0170 


- 


A++* 


1.385(51) 


1.01 


0.17 


0.9129 


0.0871 


- 


A+ + ** 


1.630(97) 


0.89 


0.07 


0.4183 


0.5817 


- 


A+ + *±* 


1.755(98) 


1.07 


0.23 


0.9563 


0.0437 


- 


A+ + **** 


1.95(16) 


0.97 


0.38 


0.9647 


0.0353 


- 


At + 


0.792(14) 


0.98 


0.03 


0.9842 


- 


0.0158 


A++* 


1.405(51) 


1.02 


0.09 


0.9249 


- 


0.0751 


A+ + ** 


1.93(14) 


1.08 


0.06 


0.9718 


- 


0.0282 


At + 


1.485(82) 


0.82 


0.06 


- 


0.8321 


0.1679 


A++* 


1.94(19) 


0.99 


0.47 


- 


0.2613 


0.7387 



Table 11: Masses in the Af + channel at the parameters shown in the header. A wide range of different 
variational sets is employed and the results are compared. In particular, we want to highlight the mixing 
with multi-glueballs and bi-torelon states. 





SU(3) 


at (3 = 


5.8945 and N L = 


= 18 




R PC 


am(a) 


\c n \ 2 


x 2 


mixc 


mixg 


mixy 


Al + 


0.799(13) 


0.98 


0.54 


0.9618 


0.0107 


0.0276 


A++* 


1.350(39) 


0.92 


0.57 


0.9819 


0.0133 


0.0048 


A+ + ** 


1.546(71) 


0.87 


1.14 


0.6752 


0.3098 


0.0150 


A+ + *±* 


1.90(13) 


0.84 


0.83 


0.9917 


0.0065 


0.0019 


At + 


0.799(13) 


0.98 


0.62 


0.9904 


0.0096 


- 


A++* 


1.347(39) 


0.91 


0.67 


0.9873 


0.0127 


- 


A+ + ** 


1.543(70) 


0.87 


1.10 


0.6851 


0.3149 


- 


A+ + *±* 


1.91(13) 


0.84 


0.83 


0.9936 


0.0064 


- 


A++ 


0.799(13) 


0.98 


0.52 


0.9801 


- 


0.0199 


A++* 


1.353(39) 


0.91 


0.60 


0.9967 


- 


0.0033 


A+ + ** 


1.92(13) 


0.85 


0.93 


0.9964 


- 


0.0036 


At + 


1.493(71) 


0.42 


0.35 




0.8056 


0.1944 



Table 12: Masses in the Af + channel at the parameters shown in the header. A wide range of different 
variational sets is employed and the results are compared. In particular, we want to highlight the mixing 
with multi-glueballs and bi-torelon states. 
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SU(4) 


at j3 = 


10.789 and N L = 


= 12 




R PC 


am{<j) 


K\ 2 


x 2 


mixc 


niixg 


mixT 


A++ 


0.808(15) 


0.99 


0.16 


0.9916 


0.0066 


0.0018 


A++* 


1.388(50) 


0.94 


0.06 


0.8399 


0.1232 


0.0369 


A++*+ 


1.605(93) 


0.93 


0.35 


0.5711 


0.1604 


0.2685 


A+ + *** 


1.607(91) 


0.88 


0.01 


0.5590 


0.4122 


0.0288 


At + 


0.808(15) 


0.99 


0.15 


0.9934 


0.0066 


- 


A++* 


1.400(51) 


0.94 


0.04 


0.8873 


0.1127 


- 


A++** 


1.596(87) 


0.88 


0.13 


0.6383 


0.3617 


- 


A+ + *** 


1.92(16) 


0.86 


0.32 


0.9878 


0.0122 


- 


At + 


0.810(15) 


0.99 


0.20 


0.9973 


- 


0.0027 


A++* 


1.404(53) 


0.94 


0.28 


0.9553 


- 


0.0447 


A++** 


1.614(90) 


0.93 


0.02 


0.6587 


- 


0.3413 


A+ + *** 


1.96(18) 


0.85 


0.03 


0.9979 


- 


0.0021 


A++ 


1.642(92) 


0.92 


0.91 




0.8334 


0.1666 



Table 13: Masses in the A{ channel at the parameters shown in the header. A wide range of different 
variational sets is employed and the results are compared. In particular, we want to highlight the mixing 
with multi-glueballs and bi-torelon states. 





SU(5) 


at j3 = 


17.107 and N L = 


12 




R PC 


am(a) 


\c n \ 2 


x 2 


mixc 


rnixg 


mixT 


A++ 


0.778(14) 


0.95 


0.01 


0.9914 


0.0059 


0.0027 


A++* 


1.409(56) 


0.93 


0.36 


0.9643 


0.0267 


0.0089 


A++** 


1.565(80) 


0.86 


0.41 


0.3586 


0.5395 


0.1019 


A+ + *** 


1.606(81) 


0.88 


0.94 


0.5318 


0.2186 


0.2496 


A+ + ***± 


1.97(16) 


0.99 


0.30 


0.9668 


0.0166 


0.0166 


At + 


0.778(14) 


0.95 


0.01 


0.9940 


0.0060 


- 


A++* 


1.402(55) 


0.92 


0.22 


0.9752 


0.0248 


- 


A+ + ** 


1.525(81) 


0.83 


0.22 


0.3203 


0.6797 


- 


A+ + *** 


1.92(15) 


0.97 


0.01 


0.9830 


0.0170 


- 


At + 


0.778(14) 


0.95 


0.01 


0.9973 


- 


0.0027 


A++* 


1.400(53) 


0.92 


0.40 


0.9876 


- 


0.0124 


A+ + ** 


1.637(85) 


0.91 


1.27 


0.6247 


- 


0.3753 


A+ + *** 


1.95(15) 


0.97 


0.10 


0.9817 


- 


0.0183 


At + 


1.565(84) 


0.87 


0.24 




0.9558 


0.0442 



Table 14: Masses in the Af + channel at the parameters shown in the header. A wide range of different 
variational sets is employed and the results are compared. In particular, we want to highlight the mixing 
with multi-glueballs and bi-torelon states. 
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SU(6) 


at j3 = 


24.845 and N L = 


= 12 




R PC 


am(a) 


K\ 2 


x 2 


niixG 


mixs 


mixT 


A++ 


0.786(14) 


0.97 


0.51 


0.9935 


0.0029 


0.0036 


A++* 


1.397(48) 


0.94 


0.57 


0.9529 


0.0367 


0.0103 


A+ + ** 


1.632(84) 


0.91 


0.39 


0.6215 


0.3170 


0.0615 


A+ + * + * 


1.650(95) 


0.90 


0.46 


0.7162 


0.0988 


0.1850 


At + 


0.786(14) 


0.97 


0.51 


0.9971 


0.0029 


- 


A++* 


1.408(49) 


0.94 


0.26 


0.9653 


0.0347 


- 


A+ + ** 


1.579(87) 


0.86 


0.02 


0.6877 


0.3123 


- 


Af + 


0.789(14) 


0.97 


0.51 


0.9963 


- 


0.0037 


A++* 


1.392(50) 


0.92 


0.41 


0.9886 


- 


0.0114 


A+ + ** 


1.71(11) 


0.96 


1.43 


0.6609 


- 


0.3391 


A++ 


1.615(89) 


0.90 


0.02 




0.9608 


0.0392 



Table 15: Masses in the Af + channel at the parameters shown in the header. A wide range of different 
variational sets is employed and the results are compared. In particular, we want to highlight the mixing 
with multi-glueballs and bi-torelon states. 





SU(7) 


at (3 = 


33.995 and N L = 


= 12 




R PC 


am(a) 


\Cn\ 2 


x 2 


mix G 


mixs 


mix<r 


A++ 


0.771(15) 


0.95 


0.33 


0.9922 


0.0048 


0.0029 


A++* 


1.553(63) 


1.06 


0.13 


0.9525 


0.0423 


0.0052 


A+ + ** 


1.658(96) 


0.94 


0.62 


0.5173 


0.4752 


0.0074 


,,++-*** 


1.82(11) 


1.01 


0.01 


0.6839 


0.0699 


0.2461 


A+ + **** 


2.05(17) 


1.09 


0.01 


0.9720 


0.0229 


0.0052 


At + 


0.771(14) 


0.95 


0.36 


0.9951 


0.0049 


- 


A++* 


1.552(63) 


1.06 


0.08 


0.9589 


0.0411 


- 


A+ + ** 


1.651(95) 


0.93 


0.29 


0.4934 


0.5066 


- 


/!+ + *** 


2.02(17) 


1.06 


0.12 


0.9776 


0.0224 


- 


At + 


0.771(15) 


0.95 


0.31 


0.9970 


- 


0.0030 


A++* 


1.562(63) 


1.06 


0.18 


0.9939 


- 


0.0061 


A+ + ** 


1.84(11) 


1.03 


0.19 


0.7515 


- 


0.2485 


A+ + **± 


2.05(16) 


1.10 


0.12 


0.9928 


- 


0.0072 


At + 


1.602(87) 


0.87 


0.06 


- 


0.8816 


0.1184 


A++* 


2.30(36) 


0.84 


0.15 


- 


0.5081 


0.4919 



Table 16: Masses in the A^ + channel at the parameters shown in the header. A wide range of different 
variational sets is employed and the results are compared. In particular, we want to highlight the mixing 
with multi-glueballs and bi-torelon states. 
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Masses in the Ax representation 




At + 


A++* 


A+ + (S) 


A+ + * + 


Ai + 


Ar+(S) 


SU(3) 


0.792(14) 


1.370(50) 


1.485(82) 


1.86(13) 


1.370(56) 


- 


SU(4) 


0.808(15) 


1.388(50) 


1.642(92) 


1.92(16) 


1.377(53 


2.54(52) 


SU(5) 


0.778(14) 


1.409(56) 


1.565(84) 


1.97(16) 


1.400(60) 


- 


SU(6) 


0.786(14) 


1.397(48) 


1.615(89) 


- 


1.315(48) 


2.43(40) 


SU(7) 


0.771(15) 


1.553(63) 


1.602(87) 


2.05(17) 


1.336(51) 




SU(8) 


0.785(16) 


1.408(63) 


- 


- 


- 


- 


SU(oo) 


0.778(8) 


1.456(41) 


1.578(47) 


2.061(28) 


1.407(17) 


2.33 


SU(oo) 


0.787(5) 


1.412(24) 


1.604(16) 


1.937(41) 


1.381(9) 


2.47(6) 



Table 17: Values of the masses in units of the lattice spacing a for each SU(JV) gauge group in the A\ 
representation. The SU(oo) masses conies from a constant fit (lower) and a constant plus 1/N 2 corrections 
(upper). Here we quote only best \ 2 fits done without using boldface values in the table. 



Masses in the A 


2 representation 




A++(2T) 


A + 2 - 


K~ 


SU(3) 


1.81(12) 


2.33(38) 


2.41(42) 


SU(4) 


1.83(15) 


2.08(19) 


2.34(37) 


SU(5) 


1.88(16) 


2.19(22) 


- 


SU(6) 


- 


2.59(51) 


- 


SU(7) 


1.84(13) 


2.43(43) 


- 


SU(8) 


- 


2.21(27) 


- 


SU(oo) 


- 


2.61(31) 


2.25 


SU(oo) 


- 


2.16(9) 


2.37(4) 



Table 18: Values of the masses in units of the lattice spacing a for each SU(JV) gauge group in the Ai 
representation. Results are obtained using only single-glueballs and bi-torelons operators. The SU(oo) 
masses comes from a constant fit (lower) and a constant plus 1/N 2 corrections (upper). Here we quote only 
best \ 2 fits done without using boldface values in the table. 







Masses in the E 


1 rcprcscnta 


don 








E++ 


E++* 


E++(2T) 


E+- 


E-+ 


E~+* 


E— 


SU(3) 


1.282(38) 


- 


1.601(82) 


2.57(44) 


1.657(95) 


2.43(38) 


- 


SU(4) 


1.290(40) 


- 


1.617(79) 


2.85(70) 


1.683(93) 


- 


2.25(25) 


SU(5) 


1.315(44) 


1.91(14) 


1.591(80) 


2.85(91) 


1.73(10) 


2.48(43) 


2.31(30) 


SU(6) 


1.293(43) 


1.82(12) 


1.620(87) 


2.80(82) 


1.78(11) 


2.56(45) 


2.20(24) 


SU(7) 


1.279(44) 


1.99(15) 


1.82(11) 


2.58(62) 


1.78(11) 


2.59(42) 


2.16(28) 


SU(8) 


1.339(53) 


1.80(14) 


1.78(12) 


2.67(72) 


1.90(14) 


- 


2.50(41) 


SU(oo) 


1.310(15) 


1.82(16) 


- 


2.714(86) 


1.830(35) 


2.589(37) 


2.35(15) 


SU(oo) 


1.296(8) 


1.871(42) 


- 


2.64(5) 


1.738(32) 


2.569(37) 


2.777(52) 



Table 19: Values of the masses in units of the lattice spacing a for each SU(iV) gauge group in the E 
representation. The SU(oo) masses comes from a constant fit (lower) and a constant plus 1/N 2 corrections 
(upper). Here we quote only best \ 2 fits done without using boldface values in the table. 
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Masses in the T\ representation 




T++ 


T+- 


T+-* 


T - + 


Tf" (2T) 


SU(3) 


2.05(21) 


1.611(78) 


- 


- 


2.35(33) 


SU(4) 


2.13(22) 


1.639(77) 


2.15(19) 


2.37(36) 


2.45(39) 


SU(5) 


2.09(21) 


1.625(76) 


1.99(16) 


2.30(30) 


2.16(28) 


SU(6) 


2.30(28) 


1.658(77) 


1.94(16) 


2.43(43) 


- 


SU(7) 


2.42(31) 


1.617(82) 


1.89(15) 


2.48(46) 


- 


SU(8) 


2.25(33) 


1.693(97) 


1.98(17) 


2.57(57) 


2.16(27) 


SU(oo) 


2.310(80) 


1.659(19) 


1.840(55) 


2.50(14) 




SU(oo) 


2.165(54) 


1.638(11) 


1.980(40) 


2.378(48) 


- 



Table 20: Values of the masses in units of the lattice spacing a for each SU(iV) gauge group in the Ti 
representation. Results are obtained using only single-glueballs and bi-torelons operators. The SU(oo) 
masses comes from a constant fit (lower) and a constant plus 1/N 2 corrections (upper). Here we quote only 
best x 2 fits done without using boldface values in the table. 



Masses in the Ti representation 




1 2 


1 2 


T +- 

1 2 


T+"(2T) 


1 2 


T 2 ~ 


SU(3) 


1.292(41) 


1.79(11) 


- 


1.87(14) 


1.665(93) 


2.26(28) 


SU(4) 


1.333(47) 


1.76(11) 


- 


1.96(16) 


1.678(95) 


2.21(25) 


SU(5) 


1.391(48) 


1.91(14) 


1.90(15) 


1.94(15) 


1.71(10) 


2.01(21) 


SU(6) 


1.263(39) 


1.99(14) 


2.08(18) 


- 


1.85(12) 


2.31(30) 


SU(7) 


1.367(50) 


1.96(15) 


2.07(17) 


2.23(23) 


1.77(11) 


2.45(33) 


SU(8) 


1.393(57) 


1.90(15) 


2.07(21) 


2.26(28) 


1.79(14) 


2.20(32) 


SU(oo) 


1.354(42) 


1.983(71) 


2.05 


- 


1.812(38) 


2.327(91) 


SU(oo) 


1.329(22) 


1.857(39) 


2.076(5) 


- 


1.73(3) 


2.299(51) 



Table 21: Values of the masses in units of the lattice spacing a for each SU(iV) gauge group in the T2 
representation. Results are obtained using only single-glueballs and bi-torelons operators. The SU(oo) 
masses comes from a constant fit (lower) and a constant plus 1/N 2 corrections (upper). Here we quote only 
best x 2 fits done without using boldface values in the table. 
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SU(oo) 


R PC 


am(a) 


c 


range N 


x 2 


At + 


0.778(8) 


0.18(0.15) 


(3,4,5,6,7,8) 


0.66 


A++* 


1.456(41) 


-0.87(0.68) 


(3,4,5,6,7,8) 


1.08 


/!+ + ** 


2.061(28) 


-1.9(0.4) 


(3,4,5,7) 


0.03 


At + (S) 


1.578(47) 


0.71(1.1) 


(4,5,6,7) 


0.17 




1.643(40) 


-1.2(0.6) 


(3,4,5,6,7) 


0.30 


A-, + 


1.407(17) 


-0.36(0.22) 


(3,4,5) 


0.04 




1.331(28) 


0.47(0.47) 


(3,4,5,6,7) 


0.41 


A r + (s) 


2.33 


3 


(4,6) 


- 


At~ 


2.61(31) 


-9(6) 


(4,5,6) 


0.29 




2.26(14) 


-1.1(2.5) 


(3,4,5,6,7,8) 


0.34 


Ar 


2.25 


1.4 


(3,4) 


- 


E++ 


1.310(15) 


-0.26(0.24) 


(3,4,5,6,7,8) 


0.23 


E++* 


1.82(16) 


2(6) 


(5,6,7,8) 


0.57 


E+- 


2.714(86) 


-1.2(1.2) 


(3,5,6,7,8) 


0.03 




2.742(99) 


-1.2(1.4) 


(3,4,5,6,7,8) 


0.04 


E~+ 


1.830(35) 


-1.8(0.6) 


(3,4,5,6,7,8) 


0.20 


E-+* 


2.589(37) 


-1.5(0.6) 


(3,5,6,7) 


0.01 


E+- 


2.714(86) 


-1.2(1.2) 


(3,5,6,7,8) 


0.03 




2.742(99) 


-1.2(1.4) 


(3,4,5,6,7,8) 


0.04 


E— 


2.35(15) 


-2(4) 


(4,5,6,8) 


0.18 




2.26(13) 


-0.2(3.2) 


(4,5,6,7,8) 


0.18 


T++ 


2.310(80) 


-2.5(1.2) 


(3,4,5,6,7,8) 


0.17 


T+- 


1.659(19) 


-0.4(0.3) 


(3,4,5,6,7,8) 


0.11 




1.621(51) 


0.1(0.8) 


(3,4,5,6,7,8) 


0.79 


T+-* 


1.840(55) 


4(2) 


(4,5,6,7,8) 


0.11 


7T+ 


2.50(14) 


-3(3) 


(4,5,6,8) 


0.06 




2.524(97) 


-3(2) 


(4,5,6,7,8) 


0.04 


1 2 


1.354(42) 


-0.5(0.7) 


(3,4,5,6,7,8) 


1.58 




1.309(22) 


-0.07(0.38) 


(3,4,5,6,7,8) 


0.50 


1 2 


1.983(71) 


-2(1) 


(3,4,5,6,7) 


0.38 




1.965(54) 


-1.9(0.8) 


(3,4,5,6,7,8) 


0.30 


1 2 


2.05 


0.8 


(6,8) 


- 




2.23(8) 


-8(3) 


(5,6,7,8) 


0.09 


rr-+ 
1 2 


1.812(38) 


-1.5(0.6) 


(3,4,5,6,7,8) 


0.23 


T 2 — 


2.327(91) 


-0.6(1.5) 


(3,6,7,8) 


0.15 




2.21(13) 


-0.1(2.2) 


(3,4,5,6,7,8) 


0.39 



Table 22: Spectrum of the SU(oo) lattice gauge theory. The masses, in units of lattice spacing, are 
obtained from fits over the range of iV shown; also the fitted coefficient c of the 1/N 2 correction is shown. 
When 2 fits are shown, the second one includes masses that we defined not reliable if the range of N is 
different. 
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Figure 4: Variational calculation for SU(3) using different sets of operators on a Nl — 12 lattice. The 
unfilled symbols represent masses that cannot be reliably interpreted as pure glueballs. 
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Figure 5: Variational calculation for SU(4) using different sets of operators on a Nl — 12 lattice. The 
unfilled symbols represent masses that cannot be reliably interpreted as pure glueballs. 
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Figure 6: Variational calculation for SU(5) using different sets of operators on a Nl — 12 lattice. The 
unfilled symbols represent masses that cannot be reliably interpreted as pure glueballs. 
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Figure 7: Variational calculation for SU(6) using different sets of operators on a Nl — 12 lattice. The 
unfilled symbols represent masses that cannot be reliably interpreted as pure glueballs. 
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Figure 8: Variational calculation for SU(7) using different sets of operators on a Nl — 12 lattice. The 
unfilled symbols represent masses that cannot be reliably interpreted as pure glueballs. 
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Figure 9: Variational calculation for SU(3) using different sets of operators on a Nl — 18 lattice. The 
unfilled symbols represent masses that cannot be reliably interpreted as pure glueballs. 
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Figure 10: Extrapolation to N — > oo of the states in the A{ channel. 
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Figure 11: Extrapolation to TV — > oo of the states in the A x + channel. We denote the less reliable states 
with open symbols. 
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Figure 12: Extrapolation to N — > oo of the states in the Af channel. We denote the less reliable states 
with open symbols. 
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Figure 13: Extrapolation to TV — > oo of the states in the E ++ channel. 
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Figure 14: Extrapolation to N — > oo of the states in the E + channel. 
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Figure 15: Extrapolation to N — > oo of the states in the E + and E channels. We denote the less 
reliable states with open symbols. 



33 




Figure 16: Extrapolation to N — > oo of the states in the T 1 ++ channel and in the T x + channel. We denote 
the less reliable states with open symbols. 




Figure 17: Extrapolation to TV — > oo of the states in the T^ channel. 
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Figure 18: Extrapolation to TV — > oo of the states in the T^~ + channel. We denote the less reliable states 
with open symbols. 
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Figure 19: Extrapolation to TV — > oo of the states in the T^~ , T 2 + and T 2 channels. We denote the 
less reliable states with open symbols. 
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Figure 20: The spectrum at N = oo. The yellow boxes represent the large N extrapolation of masses 
obtained in ref. |B8l. 
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Figure 21: Chew-Frautschi plot of the glueball spectrum 
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